o 

(N 



(N 



Observation of a Dirac point in microwave experiments with a photonic crystal 

modehng graphene 
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We present measurements of transmission and reflection spectra of a microwave photonic crystal 
composed of 874 metallic cylinders arranged in a triangular lattice. The spectra show clear evidence 
of a Dirac point, a characteristic of a spectrum of relativistic massless fermions. In fact, Dirac points 
are a peculiar property of the electronic band structure of graphene, whose properties consequently 
can be described by the relativistic Dirac equation. In the vicinity of the Dirac point, the measured 
reflection spectra resemble those obtained by conductance measurements in scanning tunneling 
microscopy of graphene flakes. 

PACS numbers: 42.70.Qs, 73.22.Pr, 42.25.Fx 
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Graphene is a monolayer of carbon atoms arranged in 
a honeycomb lattice [l|, |2| . Due to its peculiar electronic 
properties this carbon allotrope recently attracted a lot of 
attention in condensed matter physics. The conduction 
and the valence band of the electronic energy in graphene 
form conically shaped valleys that touch each other at 
the corners of the Brillouin zone [3| . As a consequence, 
close to these touch points the energy of the electron de- 
pends linearly on its quasi-momentum vector. This linear 
dispersion relation implies an energy independent veloc- 
ity, and the related wave equation is the Dirac equation 
[J|. Thus, although the Fermi velocity is typically 300 
times smaller than that of light, the energy spectrum of 
the electrons in graphene is similar to that of massless 
relativistic fermions [^-iTJ. However, graphene with its 
Dirac spectrum is not an exception. Photonic crystals, 
an optical analogue of ordinary crystals, possess similar 
particular properties [8|, |9|. The unusual transmission 
properties near a Dirac point predicted in [9!| were ob- 
served experimentally both in sonic and in microwave 
photonic crystals |10l.lll|. 

The photonic crystal considered in the present work is 
two-dimensional and composed of rows of metallic cylin- 
ders which are arranged to form the triangular lattice 
schematically shown in the left part of Fig. [T] Electro- 
magnetic waves propagating in such a periodic structure 
exhibit a dispersion relation with a band structure sim- 
ilar to the electronic band structure in a solid. That 
for the triangular lattice of metallic cylinders with ra- 
dius R = 0.25a, where a is the lattice constant, is shown 
on the right part of Fig. [TJ It was obtained by solving 
numerically the Helmholtz equation with the finite dif- 
ference method [1^. In this so-called band diagram the 
frequency / is plotted as function of the quasi momentum 
components {k,j;, ky) for the first two propagating modes. 
Around the corners of the first Brillouin zone the bands 
have the shape of cones and the band structure resembles 
that of the electronic energy in graphene [2, |5[ . In fact 
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FIG. 1: (Color online) Left: Triangular lattice of metalUc 
cy Anders. The arrows indicate the two different directions FM 
and FK in the triangular lattice; the radius R of the metallic 
cylinders equals R = 0.25a, where a is the lattice constant. 
Right: The plot of the numerically determined band structure 
f{kx,ky), as function of the quasi momentum components 
(kx, ky). The solid lines indicate the hexagonal Brillouin zone. 



both systems have the same Bravais lattice and the wave 
functions with quasi momentum at the corners of the first 
Brillouin zone possess the same symmetry group [13| . In 
the hterature the touching point is referred to as Dirac 
point [2|, |lj| and its frequency as the Dirac frequency. 

In the vicinity of a touching point the Hamiltonian 
governing the spectrum can be written as a 2 x 2 matrix 
in the basis of the doublet states which are degenerate at 
the Dirac frequency and has the form [14j 



H = ujot + uc {5k^d^ 



5kydy) 



(1) 



Here, Ux, Oy are the Pauli matrices, 8k = ((5/ca:, 8ky) is the 
displacement vector with respect to a corner of the first 
Brillouin zone, w — 2nf, with / denoting the excitation 
frequency, ujd is the Dirac frequency and vd = |Vi^(^)| is 
the group velocity which is approximately constant close 
to the Dirac point. Diagonalization of H in Eq. ([1]) yields 



w(|(5fc|) = uj£) ± V]j\/Skx + 5k 



(2) 



which reproduces the hnear dispersion relation around 
the tips of the cones of the Brillouin zones shown in Fig.[T] 
The substitution Sk^ -^ 
leads to the Dirac equation 
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idx and 6ky — !• —idy in Eq. ([T]) 





dx + idy 
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(3) 

where V'l, V'2 represent the amplitudes of the degenerate 
doublet state at the Dirac frequency. Thus, despite of 
the fact that photons are bosons in the vicinity of ujd the 
waves in a photonic crystal can be effectively described 
by the Dirac equation for fermions with spin 1/2. 

The similarity of the band structure of a macroscopic 
photonic crystal with the electronic band structure of 
graphene, which is experimentally much more difficult to 
access, allows the experimental study of various relativis- 
tic phenomena. The photonic crystal used in the exper- 
iments is composed of a total of 874 metallic cylinders 
with radius i? = 5 mm and height h = 8 mm squeezed 
between two metallic plates. The experimental setup is 
shown in Fig. [2] Each cylinder is screwed to the top and 
bottom copper plate to ensure a proper electric contact 
even at high frequencies and thus reproducibility of the 
measurements. The resulting photonic crystal consists of 
38 X 23 cylinders and has a size of 400 x 900 mm. Below a 
certain excitation frequency fc — c/2h, with c being the 
speed of light and h the gap between the plates, only the 
lowest transverse magnetic (TM) mode with the electric 
field perpendicular to the plates can propagate. Con- 
sequently, like for two-dimensional microwave billiards 
[ISl |l6| , the propagating modes correspond to solutions 
of the two-dimensional scalar Helmholtz equation, which 
is mathematically identical to the Schrodinger equation 
for the corresponding quantum multiple-scattering prob- 
lem. 

To determine the band structure of the photonic 
crystal, transmission measurements are performed with 
nearly plane waves. To produce and detect them a two- 
dimensional horn antenna a and a wave-guide antenna 
b are placed, respectively, on the opposite sides of the 
crystal 240 mm apart from it (see Fig. [2]). Both are 
screwed tightly to the top and bottom plate and are at- 
tached to waveguide-to-coaxial cable adapters. Trans- 
mission (reflection) spectra are measured with an Agi- 
lent PNA-L 5230 vectorial network analyzer (VNA). Mi- 
crowave power is emitted into the region of the photonic 
crystal via one antenna a and the magnitude and phase 
of the transmitted (reflected) power at the other (same) 
antenna b (a) is determined relative to the input signal, 
thus yielding the complex valued scattering matrix ele- 
ment Sbaif)- Systematic errors due to the attenuation of 
microwave power in the coaxial cables and the reflection 
at the connectors were removed by a proper calibration 
of the VNA. Ribbons of microwave absorption foam were 
placed between the plates at their edges to avoid disturb- 
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FIG. 2: Scheme of the experimental setup. The photonic 
crystal consists of altogether 874 metallic cylinders arranged 
in a triangular lattice. Electromagnetic waves are transmitted 
from the horn antenna through the lattice to the wave guide 
antenna. The antennae and the lattice are oriented with re- 
spect to each other such that the waves impinge the lattice in 
TM direction (see Fig[T]). 



ing reflections. 

The photonic crystal is arranged such, that the waves 
produced by the horn antenna impinge it in FM direction 
(see Figs. [T] and [5|). The left part of Fig. [3] shows a mea- 
sured transmission spectrum \Sbaif)\'^ through the pho- 
tonic crystal in FM direction. It shows two stop bands, 
one below 11.5 GHz and one between 18 and 19.8 GHz, 
where the transmission is reduced by 4 orders of magni- 
tude. In between we observe transmission with a broad 
minimum around 14 GHz. The right part shows the cal- 
culated band structure along the path FMKF (see in- 
set of Fig. [S]) inside the first Brillouin zone. The first 
and second band touch at the corner K of the Brillouin 
zone. The positions of the experimental stop bands are 
in good agreement with the band gaps in the calculated 
band structure. The minimum in the transmission spec- 
trum is located around the Dirac frequency and indicates 
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FIG. 3: Left: The measured transmission spectrum from the 
horn to the wave guide antenna. Right: Numerically calu- 
lated band structure /(|fc|) with k chosen along the FMRF 
path. The inset shows the Brillouin zone of the lattice and 
the irreducible Brillouin zone FMK. 



a partial band gap [17] between the first and second band 
expected for the transmission of plane waves in FM di- 
rection. However, since the incident wave is not perfectly 
plane and thus contains components with propagation di- 
rection different from the FM one, this band gap is not as 
well pronounced as the other two. Altogether the exper- 
imental transmission spectrum resembles the predicted 
band structure of the photonic crystal. 

A very effective method to localize a Dirac point is to 
measure the frequency dependence of the local density 
of states (LDOS). It is defined as the imaginary part of 
the Green function for the photonic crystal, L{f, f) — 
— l/7rlm G{r,r,f). According to [l8| it can be written 
in terms of the frequencies fn{k) and the wave functions 
ipn{k, r) of the nth propagating mode as 

L{r,f) = -^ [ y^\i;^{k,r)\'^5{f - Uk))cPk, 

where the integration over the quasi momentum k ex- 
tends over the first Brillouin zone with area Abz- Us- 
ing the dispersion relation Eq. 1^ and that the modulus 
squared of the wave function depends only weakly on k 
close to the Dirac point yields 
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(5) 



Here, (IV'(^a)l^) denotes the average over k of the mod- 
ulus squared of the wave functions in the vicinity of the 
Dirac point. Equation ([5]) reflects the well-known theo- 
retical result [si y, [l3| that near the Dirac frequency the 
LDOS tends to zero linearly in /. The LDOS can be de- 
termined from the measured reflection spectra \Saa{f)\^ 
with the antenna placed at a position Tq within the crys- 
tal. The connection between the LDOS and the measured 
reflection spectra is established via the scattering matrix 
formalism for microwave resonators 1&, 20' originally de- 
veloped in the context of compound-nucleus reaction the- 
ory 21|. Within this formalism the scattering matrix 
can be written in the form S{f) = (1 — iK){'t + iK)~^, 
where K = W^{E — H)^^W with H denoting the Hamil- 
tonian for the photonic crystal and W the matrix cou- 
pling the open channels, as e.g., the antenna channel 
to the modes of the photonic crystal. Expressing K 
in terms of the Green function of the photonic crystal, 
G{ra,rb,E + ie)=J2n^nin)Mra){E + ie-En)~^ leads 
to the relation 



l-|5aa(/)p0C 



Lira J) 



\G{ra,raJ) + A\ 



(6) 



Here, Gfry,, tv,., f) denotes the renormalized Green func- 
tion [22!] and the constant A depends on the coupling of 
the antenna to the photonic crystal. In the vicinity of the 
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FIG. 4: Left: A reflection spectrum measured with the wire 
antenna located in the middle of the crystal in the center of a 
triangle of cylinders as schematically shown in the inset. The 
sharp resonances at the edges of the bands are related to the 
so-called van Hove singularities. Right: For camparison the 
computed band diagram from Fig. [3] is shown. 



Dirac point one can neglect the frequency dependence of 
the denominator, i.e., replace it by a constant B{fa, fo) 
and finally gets the relation between the experimental 
observable \Saa{f)\'^ and the LDOS 



\Saa{f)\^^l-B{rajD)LiraJ)- 



(7) 



Thus, according to Eq. ([5]), the reflected power \Saa{f)\^ 
tends linearly to unity in the vicinity of the Dirac fre- 
quency. This is the "fingerprint" used to identify the 
Dirac point. 

For the measurement of the reflection spectra a wire 
antenna consisting of a metallic pin of 1 mm in diameter 
was lead through a 3 mm wide drilling in the top plate. 
The left part of Fig. |4] shows the reflection spectrum mea- 
sured with an antenna placed in the middle of the pho- 
tonic crystal as shown in the inset of this figure. The 
location of the antenna is chosen in the center of three 
cylinders forming a triangle to minimize the disturbance 
of the propagating mode at the Dirac frequency. The ex- 
perimental reflection spectrum has a clearly pronounced 
maximum around «14 GHz, i.e., within the frequency 
range where the Dirac point is expected, and shows the 
characteristic cusp structure. The sharp resonances at 
the edges of the bands are related to the so-called van 
Hove singularities [23J. 

The description of the experimental reflection spec- 
trum around the Dirac frequency based on Eqs. ([5|) and 
d?]) allows to express the measured quantity \Saa{f)\'^ in 
terms of a 3-parameter formula. 



\Saa{f)\'^D-C\f-fDl 



(8) 



where the parameter D takes into account that due to 
Ohmic losses the reflection deviates from unity even in 
the case of vanishingly small LDOS. The parameter C 
describes the slope of the cusp in the reflection spectrum 
close to the Dirac frequency. The upper and the lower 




FIG. 5: Two measured reflection spectra in the vicinity of 
the Dirac frequency. The upper spectrum was measured with 
an antenna placed in the middle of the crystal, the lower one 
with an antenna placed 6 rows apart from the boundary of 
the crystal. The smooth solid lines are fits of Eq. ([S| to the 
experimental spectra. 



parts of Fig.[5]show the experimental reflection spectra in 
the frequency range 11-17 GHz for two different antennae 
together with the fit of Eq. ([8]) to the data. The spec- 
trum shown in the upper panel of Fig. [S] was measured 
with an antenna placed in the middle of the crystal, that 
shown in the lower panel was obtained with one placed 
6 rows apart from the crystal boundary. In the former 
case the antenna reached 2 mm into the space between 
the plates, in the latter only 0.5 mm. The upper spec- 
trum is well described by Eq. ([S]) and the values of the fit 
parameters are D = 0.826±0.003, C = 0.413±0.006, and 
the Dirac frequency Jd = 13.797 ± 0.004 GHz is in good 
agreement with the calculated value fo = 13.81 GHz. 
The group velocity vd at the Dirac point can not be ex- 
tracted from C because the value of (|'0(^a)P) in Eq. ([5]) 
is unknown. The shape and position of the cusp are re- 
covered in the lower reflection spectrum in Fig. [5] The 
obtained values of the parameters are D — 0.988 ±0.001, 
C = 0.329 ± 0.003 and fo = 13.793 ± 0.002 GHz. The 
value for the Dirac frequency differs slightly from that 
determined for the antenna placed in the middle of the 
crystal but the difference is within the standard error of 
the fit. The value of the parameter D is closer to unity 
in the latter case because the antenna reaches less into 
the space between the plates and thus couples weaker 
to the propagating modes. This is also reflected in the 
slope C of the cusp, which is not as steep. Moreover, the 
cusp at the Dirac frequency is rounded in the lower spec- 
trum. This is attributed to the proximity of the antenna 
to the crystal boundary. The fluctuations observed in 
the reflection spectra around the mean behaviour given 
by Eq. ([5]) are attributed to quasi bound modes trapped 
in the photonic crystal. In the vicinity of the Dirac point, 
the measured reflection spectra closely resemble those ob- 



tained by conductance measurements in scanning tunnel- 
ing microscopy of graphene flakes deposited on a graphite 
substrate (24 |. 

In summary, we have measured the transmission and 
the reflection spectra of a photonic crystal with a tri- 
angular lattice composed of metallic cylinders. A cusp 
structure was observed in the reflection spectra close to 
the expected Dirac frequency and related to the local 
density of states in the photonic crystal, providing clear 
evidence for the existence of a Dirac point. The exper- 
imental setup presented in this work is well suited for 
the study of effects connected with so-called edge states 
in graphene IJ, |23l and the predicted pseudo-diffusive 
transmission at the Dirac point. However, the most di- 
rect application is the experimental investigation of prop- 
erties of the eigenvalues and eigenfunctions of a Dirac 
billiard, by inserting the photonic crystal setup into a 
closed resonator box [26|. 
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